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Abstract 

Catalan numbers arise in many enumerative contexts as the counting 
sequence of combinatorial structures. In this work, we consider natural 
Markov chains on some of the realizations of the Catalan sequence. While 
our main result is in deriving an 0(n 2 log n) bound on the mixing time 
in L 2 (and hence total variation) distance for the random transposition 
chain on Dyck paths, we raise several open questions, including the opti¬ 
mality of the above bound. The novelty in our proof is in establishing a 
certain negative correlation property among random bases of lattice path 
matroids, including the so-called Catalan matroid which can be defined 
using Dyck paths. 


1 Introduction 

There are several combinatorial structures whose enumeration is given by the 
Catalan sequence: 1,2,5,14,42,..., where the nth term (for n > 1) is given 
by the Catalan number C n = ( 2 . Some examples include the set of all 

triangulations of a regular polygon of n + 2 sides, the set of non-crossing par¬ 
titions of an n-set (the lattice on which is of much interest to researchers in 
free probability), and the set of balanced parentheses, with n left and n right 
parentheses—or equivalently, described as Catalan strings x £ {±l} 2n of n l’s 
and n —l’s with non-negative partial sums, '£2i=i x i — for a U 3- This last 
structure (of Catalan strings) is also known as the set of Dyck paths, and is 
visualized as lattice paths of n up-steps and n down-steps, representing the l’s 
and the —l’s respectively. 

For any given n > 1, generating a Catalan structure (such as triangulation 
or a Dyck path) uniformly at random from the set of C n many, is a straight¬ 
forward task, and can be done in time linear in n, not too unlike generating a 
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permutation uniformly at random (out of the n! many permutations) of n dis¬ 
tinct letters. However, the study (of convergence to equilibrium) of particular 
Markov chain Monte Carlo algorithms which yield a random Catalan structure 
is seemingly much more interesting. The inspiration to this endeavor in part 
stems from an open problem of David Aldous (T), who conjectured that the ran¬ 
dom walk on triangulations of a polygon on n sides, performed using uniform 
diagonal flips, ought to take time roughly n 3 / 2 (up to factors logarithmic in n). 
Despite much effort by various researchers in this topic over several years, the 
best known bounds for the mixing time of the chain on triangulations remain 
those of McShine-Tetali m , who proved an upper bound of 0(n 4 ) on the re¬ 
laxation time (also known as the inverse spectral gap), and Molloy-Reed-Steiger 
m who showed fl(n 3 / 2 ) as a lower bound. 

In the work, we explore a (different) Markov chain on a different realization of 
the Catalan number: we consider the set of Dyck paths, described using binary 
strings x G {±l} 2n , as mentioned above, and consider the natural local move of 
exchanging aq and Xj, for a pair (i,j) chosen uniformly at random from among 
all possible pairs, provided that the resulting string is a (valid) Catalan string; 
else we reject such a move. Note that, without the rejection step, we would have 
a random transposition chain on the set of (many binary strings with an 
equal number of l’s and — l’s. It is well-known (see e.g., [7])—sometimes under 
the name Bernoulli-Laplace model —that 0(nlog?r) such random transpositions 
are necessary and sufficient to reach (close to) equilibrium. It is also well-known 
that the so-called random transposition shuffle of n distinct cards mixes in 
O(nlogn) time. However, the “sea-level” constraint of having the partial sums 
of the (Catalan) string stay non-negative seems to thwart any type of straight 
forward (or otherwise) analyses, despite the best efforts of several experts in the 
field. 

Our main result is an 0(n 2 log n) upper bound on the mixing time of the 
random transposition walk on Dyck paths of length n. The proof relies on a 
rephrasing of the walk as a basis exchange walk on a balanced matroid. In 
Section [2] we describe the Catalan matroid (due to Ardila [5]) and the prior 
work mmm regarding mixing times for the basis exchange walk on balanced 
matroids. All that remains to bound the mixing time of our Dyck path walk is 
the proof in Section [3] that the Catalan matroid is balanced. In Section [4] we 
give a short proof of a mixing time lower bound for Dyck random transpositions 
and, finally, in Section [5] we explore the possible consequences of this result for 
other walks on Catalan structures and pose some open questions. 

1.1 Random walks on lattice paths 

Throughout, for integers a < b we let [a, b] := {a, a + 1,... 6} denote the discrete 
closed interval and [a] := [l,a]. 

A lattice path of length m is a string P G {±l} m of up-steps (+) and down- 
steps (— ). The height of P at index i is hi(P ) = Y^j=i Pi; and we can draw P 
on the grid as the graph of fp(i) = hi(P). That is, we draw a path starting 
at (0,0) and taking up-steps (1,1) and down-steps (1,-1). We will refer to a 
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lattice path P of length m as a lattice path from (0,0) to ( m , h m {P)). 

From this we can define a partial order on the set of lattice paths from (0, 0) 
to (to, 2r — to) by letting P < Q whenever hi(P ) < hi{Q) for all i £ [m\. Note 
that P < Q if and only if qi < pi for all i £ [r], where pi (resp. qi) is the index 
of the *th up-step in P (resp. Q). 

Definition 1. A Dyck path of length 2 n is a lattice path P from (0,0) to (2n, 0) 
with hi(P ) > 0 for all i £ [2n]. 

That is, the Dyck paths are precisely the paths P from (0, 0) to (0,2 n) with 
P > (-1—)". It is a standard result that the number of Dyck paths of length 2 n 
is the Catalan number C n = ( 2 "). 

A natural class of random walks on lattice paths from (0, 0) to (to, h) is 
the transposition walk , which at each step picks random indices i,j £ [m] and 
swaps the steps of P at those indices. Indeed, if i and j are chosen uniformly 
independently at random this is just the Bernoulli-Laplace model on , where 
r = (h+m)/ 2 is the number of up-steps in every such path, and has been studied 
extensively. If we pick i. j to be uniform among adjacent pairs of indices, we 
get a different walk, whose mixing properties were determined in [16] . 

But what if we consider lattice paths from (0, 0) to (2 n, 0) and restrict the 
walk to Dyck paths? That is, from any Dyck path we pick i, j according to one 
of these rules but only perform the move if the resulting path is again a Dyck 
path. (We will call these walks the Dyck random transpositions walk and Dyck 
adjacent transpositions walk , respectively.) 

Indeed, the analysis of !l6 extends also to this case, yielding the same up¬ 
per bound for the mixing time of the Dyck adjacent transpositions walk. On 
the other hand, the Dyck random transpositions walk has evaded such precise 
analysis. 


2 Lattice path matroids 

2.1 Matroids and the basis exchange walk 

Recall that a nonempty set B C 2 U is the set of bases of a matroid M = (C7, B) 
if the following basis exchange axiom holds: 

Matroid Basis Exchange Axiom. For any bases A,B £ B and 

every e £ A\B there exists f £ B \ A such that A \ {e} U {/} € B 
is a basis. 

Among other things, this axiom guarantees that all bases have the same cardi¬ 
nality, which is the rank of M. 

We will also make use of two dual operations on matroids: contraction and 
deletion. 

Definition 2. For a matroid M = (U,B) and an element e £ M, the matroid 
M contract e is M e = (U, B e ), where B e = {B £ B : e £ B}. Similarly, M delete 
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e is M e = (U, B e ), where B e = {B £ B : e ^ B}0 

A minor of M is any matroid which can be obtained from M through a 
series of contractions and deletions. 

The order in which contractions and deletions are performed does not mat¬ 
ter, so we will write Mj for the matroid obtained from M by contracting the 
elements in I and deleting the elements in J, and £>/ will denote the set of bases 
of M{. 

Given a matroid M , Feder and Mihail [5] study the following basis exchange 
walk on the state space B of bases: 

From state B £ B, pick uniformly and independently at random 
elements a £ U and b £ B and move in the next step to B' = 

B U {a} \ {6} if B' £ B (else remain at B). 

The matroid basis exchange axiom guarantees that this walk is ergodic, and 
since it is symmetric its stationary distribution is uniform. It is tentatively 
conjectured that the basis exchange walk is fast for any matroid (i.e., the mixing 
time is bounded by some polynomial in m = |t/|), but there is little evidence in 
favor of this. 

On the other hand, [8] introduced the notion of balanced matroids, which 
capture the notion that for a randomly chosen basis (of the matroid or any of its 
minors), conditioning on the occurrence of one element (in the basis) makes the 
occurrence of any other less probable. They show that for the case of balanced 
matroids the walk is indeed rapidly mixing, by using decomposition techniques 
to bound its spectral gap. 

Definition 3. A matroid M = ( U,B) is negatively correlated if for every pair 
of distinct elements e, / £ U 


\B\ - \Bf\ • 1 j 

The matroid M is balanced if M and all of its minors are negatively correlated. 

Negative correlation is equivalent to the very natural condition that for a 
uniform random basis B £ B, Pr[e £ B] > Pr[e £ B\f £ B\. 

Indeed, many common classes of matroids are balanced, including uniform 
matroids (whose bases are all size-r subsets of the ground set), graphic matroids 
(with the ground set being the edges of a connected graph and the bases be¬ 
ing spanning trees of the graph), matroids of rank < 3, and regular matroids 
(matroids which can be represented over every field) [U [5] . We will rely on the 
following extension of Feder and Mihail’s result for balanced matroids, due to 
Jerrum and Son. 

1 Note that this definition differs slightly from the usual one, in which the element e being 
contracted or deleted is removed from the ground set. In our case it will be convenient to leave 
e in place, so as to more easily identify B e C B and preserve identities such &s B = B e U B e 
(with Be r\B e = 9). 
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Figure 1: The Dyck path (+ + + H-1-M-1-) above corresponds 

to the basis {1, 2, 3,4, 6,8, 9,13} of the Catalan matroid of order 8. 


Theorem 4 (0). The spectral gap A and log-Sobolev constant a for the basis 
exchange walk on a balanced matroid M of rank r on a ground set of size m are 
lower bounded by 

, 2 , 1 

A >- and a > — -. 

mr 2 mr 

In particular, via standard mixing time bounds (see [5]), the bound on the 
log-Sobolev constant implies that the mixing time of the basis exchange walk 
on a balanced matroid is at most 0(mr log log |£>|). 

Remark 5. Due to these standard bounds and the result for balanced matroids 
the proof of our main result never explicitly refers to the definitions of A and 
a. These definitions can, however, be found in Section U where they are used 
to give a lower bound on the mixing time for the Dyck random transposition 
chain. 

2.2 The Catalan matroid and other lattice path matroids 

To aid our setting of the Catalan transposition walk in terms of matroids, we 
refer to an observation of Ardila [2] that the set of Dyck paths can be thought 
of as a matroid. 

Definition 6. The Catalan matroid of order n is C{n) = ([2 n\,B(n)), where 
the elements of B(n) are the index-sets of up-steps in Dyck paths of length 2 n. 

Ardila notes that this is precisely the transversal matroid for the set system 
S = {[1], [3],..., [2 n — 1]} (i.e., its bases are the systems of distinct represen¬ 
tatives of S). In short, the representative of the set [2* — 1] will be the index 
of the ith up-step in the corresponding Dyck path (although this assignment of 
distinct representatives need not be unique). Indeed, 

Observation 7. The basis exchange walk on the Catalan matroid is exactly the 
random transposition walk on Dyck paths. 

In other words, given Theorem^ to obtain a mixing time bound for the Dyck 
transposition walk it suffices to show that the Catalan matroids are balanced. 

To this end, it would be convenient if the Catalan matroids belonged to some 
known class of balanced matroids. As noted above, the main class of matroids 
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known to be balanced is regular matroids, but Ardila notes that C(n) is not 
representable over any F g for q < n — 2 and thus is not regular. Transversal 
matroids also need not be balanced in general: Choe and Wagner m give a 
transversal matroid of rank 4 which is not balanced. 

Thus the bulk of our work here will be to show that the Catalan matroid is 
balanced, and from this the main mixing result will follow immediately. For our 
later discussion it will help to work with a minor-closed class of matroids (which 
the Catalan matroids certainly are not). Bonin and de Mier ([3]) discuss the 
following class of lattice path matroids, which generalize the Catalan matroids 
by allowing any pair of bounding paths. 

Definition 8 . For two lattice paths A < B from (0, 0) to (m, 2r — to), consider 
the set L of lattice paths P from (0, 0) to (to, 2 r — to) with A < P < B. The 
lattice path matroid C[A,B\ (which is of rank r on ground set [to]) has as its 
bases the index sets of up-steps of paths in L. 

Although it is not immediately obvious that C[A, B] is a matroid, [3] ob¬ 
serves that in fact the lattice path matroid C[A, B] is the transversal matroid 
of the set system {[ai, &i], [ 02 , £> 2 ], • ■ ■, \a r , &r]}, where ai (resp. 6 j) is the in¬ 
dex of the ith up-step in A (resp. B). Indeed, it is shown in that paper 
that the transversal matroids of a set systems {[ai, t>i], [< 22 , £> 2 ], ■ ■ •, [o r , b r ]} with 
a 1 <02 < • • • < a r and &i < 62 < • • • < b r are precisely the lattice path 
matroids. (The set system of this form corresponding to a given lattice path 
matroid may not be unique; for example, the Catalan matroid is the transver¬ 
sal matroid for both {[1], [3],..., [2n — 1]} and {[1], [2,3],..., [n, 2 n — 1]}.) The 
class of basis exchange walks on lattice path matroids also includes the (un¬ 
constrained) Bernoulli-Laplace model as the basis exchange walk on a uniform 
matroid. 

In addition to showing this correspondence, Bonin and de Mier show that 
the class of lattice path matroids is closed under taking minors and duals. They 
also define a smaller (in fact, minimal) minor-closed class of matroids containing 
C n , which they call generalized Catalan matroids , consisting of the lattice path 
matroids M(A,B) where A = (+) r (—) m_r is maximal among all paths from 
(0,0) to (to, 27' — to). They give several nice properties of generalized Catalan 
matroids, which we will not go into here as our analysis covers all lattice path 
matroids. Indeed, Sohoni Til has actually already shown the balanced property 
of generalized Catalan matroids, although he calls them Schubert matroids and 
does not link his result to rapid mixing of the Dyck random transposition chain 
(although he does mention the Dyck adjacent transposition chain). Our proof 
applies to the more general class of lattice path matroids. 

Given that the class of lattice path matroids is minor-closed, it suffices for 
our main result to show that lattice path matroids are negatively correlated, 
but it is worth noting that the same method can also be used to show negative 
correlation directly for any minor of a lattice path matroid. 
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3 Mixing bound for random transpositions 

To show negative correlation for lattice path matroids, it will be convenient 
to use an equivalent formulation of 0, which is easily obtained by repeatedly 
applying the identity \B\ = \B e \ + \B e \, 

\B ef \\B e f\<\Bl\\B}\. (2) 


Now we are ready to prove 

Theorem 9. For every pair of lattice paths A < B from (0, 0) to (to, 2r — m) 
the lattice path matroid C[A , B] is negatively correlated. 

Proof. Let £[A, B] = (jm],B) be any lattice path matroid. 

To prove inequality 0 for every pair e < f £ [to], we will construct an 
injective map 

<p e f : B ef x B Kf -> B{ x B). 

Note that we can associate members of these sets with lattice paths in C[A 1 B] 
by 

• B e f. paths with up-steps at indices e,/, 

• B e *: paths with down-steps at indices e,/, 

• B[: paths with an up-step at index e and a down-step at index /, and 

• By. paths with a down-step at index e and an up-step at index /. 

Let P £ B e f and Q £ B e f be lattice paths and consider the following cases. 
(The figures shown are for the Catalan matroid.) 

Case 1 . Suppose the paths P, Q intersect (without necessarily crossing) in the 
region (I) between e and /. Note that this includes all cases where the 
path P is below Q at e and above Q at / or vice versa. Take the first such 
intersection point x in (I) and switch the paths P, Q after x to obtain new 
paths P' £ B[, Q' £ B*j as shown in Figure [ 2 ] Set <p e f(P, Q ) = (P / , Q')- 



Figure 2: The injection for Case 1. 

Case 2 . Suppose the paths do not meet in region (I) and consider the paths 
P, Q in the region (II) after position /. Imagine translating the fragment 
of path Q after / to P so that their initial points in / coincide. If the 


7 


imaginary fragment intersects P, let x be the first such point of intersec¬ 
tion. Construct new paths P', Q' by swapping the segments of P and Q 
between / and x (including / itself). See Figure[3] If the new paths P', Q' 
are both between A and B (so that P' £ and Q' £ Bj), then we set 

PefiP, Q) = {P 1 , Q')- 



Figure 3: The injection for Case 2. 


Case 3. Otherwise, perform the mirror image of the operation from Case 2 
in the part (III) up to position e, as shown in Figure 0] and again set 
<p ef (P,Q) = {Q',P'). 




Figure 4: The injection for Case 3. Case 2 fails because the translated segment 
of Q leaves the valid region before intersecting Q. 


Case 1 covers all pairs of paths P, Q such that P is above Q at e and below 
it at / or vice versa. Case 2 covers all remaining pairs of paths with P above Q 
at /, since in this case the fragment of Q must intersect P and both resulting 
paths ( P',Q') are between P and Q. Finally, Case 3 covers all still remaining 
pairs of paths with P below Q at e. Hence every pair of paths P, Q is covered 
by one of these cases. 

Finally we must argue that <p e f is injective. Given P 7 , Q' we must be able 
to determine which case above was applied to produce them, and from this it 
is simple to recover P and Q. Since we know e, / it is easy to identify the 
regions (I), (II), (III) in the three cases. If P',Q' intersect in region (I) then 
they must have come from Case 1, as neither of the later cases can produce such 
an intersection. If they do not, first try to apply the inverse operation for Case 
2 , which is the same as the forwards operation; if this was not possible for the 
starting paths P, Q then it is not possible for P', Q' either (as performing the 
transformation of Case 3 cannot cause the transformation of Case 2 to become 
valid if it was not already), so if we have a pair of paths P ', Q' for which the 
second move is possible and does result in lattice paths between A and B then 


we know we must have arrived at it through Case 2. Otherwise, we were in Case 
3. See Figure [5] 



Figure 5: An example of recovering the case from the result. Since the paths 
do not intersect in region (I) we must have come from Case 2 or 3. Case 2 fails 
because the translated segment of P' leaves the valid region before rejoining Q ', 
so the paths must have come from Case 3. Note that the transformation <j> e j is 
usually not onto. 


□ 

Remark 10. Since Case 1 can be applied whenever P is above Q at e and below 
at / or vice versa, Case 2 can be applied whenever P is above Q at both e 
and /, and Case 3 can be applied whenever P is below Q at both e and /, it 
is tempting to try to divide the cases more naturally according to the relative 
heights of P and Q at e and /. However, note that in Cases 2 and 3 (and 
sometimes Case 1) the resulting paths P', Q' always have P' above Q' at both 
e and /, and this more natural division of cases does not allow us to recover 
uniquely which transformation was applied. 

Remark 11. While the result of [3] that lattice path matroids are minor-closed 
shows that it suffices for us to prove negative correlation, it is worth noting that 
the same injection shows negative correlation directly for any minor of a lattice 
path matroid. In particular, if we assume that P and Q above are both in the 
minor Bj (i.e., have up-steps at the indices of I and down-steps at the indices 
in J) then the resulting paths P' and Q' again have this property. Moreover, 
the injection also implies that lattice path matroids are Rayleigh matroids (a 
class of matroids studied in 0)- 


4 Lower bound 


To derive a lower bound on the mixing time of the Dyck random transposition 
chain we will need to refer to the definitions of the spectral gap and log-Sobolev 
constant. Recall that for a Markov chain with state space D, transition matrix 
P, and stationary distribution n 


A := inf 
/ 


£(/> /) 

Var n (f) 


and 


:= inf 
/ 


MJ) 

CM) ’ 
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where 


£(/,/) = \ - f(y)) 2 p ( x ,y)H x ) and 

x,y 

Ar(/) = 5D /( a; ) 2 * ( lo e/C 37 ) 2 - lo gE(/ 2 ))7r(a:) 

- 7T 

X 

and the infima are taken over non-constant functions / : fl —> R. 

Using standard results we can (asymptotically) lower bound the mixing time 
by the relaxation time . 

Theorem 12. The spectral gap A and log-Sobolev constant a for the Dyck ran¬ 
dom transposition chain satisfy 

4 

2a < A < 

n 

Proof. It is always the case that a < A/2 (see, e.g., [SJ), so it suffices to show 
that 1/A > n/4. 

For x G Q a Dyck path of length 2 n, consider the function f(x) giving the 
number of down-steps of x at even indices. The number of Dyck paths x of length 
2 n with f(x) = k is precisely the Narayana number N(n,k) = d ()!) (see 
[15]). Note that / is 1-Lipschitz, in that if x, y differ by a single transposition (i.e. 
P(x,y ) > 0) then | f(x) — f(y) \ < 1. We can therefore bound A by noting that 
for a 1-Lipschitz function we have £(/, /) < 1/2, and so it suffices to give a lower 
bound on Vai•„■(/), which is precisely the variance of the Narayana distribution 
with p.m.f. p(k) = N(n,k)/C n . This distribution is liypergeometric, and its 
variance is > n/8 (see, e.g., [H]). □ 

Corollary 13. The (total variation) mixing time of the Dyck random transpo¬ 
sition walk is at least D(n). 


5 Open questions 

Of course, the most immediate open question is whether the 0(rr logn) mixing 
time bound for Dyck transpositions is tight. The general result for lattice path 
matroids is clearly not tight in every case (as it gives only order n 2 log n mixing 
time for Bernoulli-Laplace, where we know that the true bound is order n log n), 
but there are lattice path matroids for which the 0(n 2 log n) bound is tigh1@. We 
believe that the Catalan matroid falls closer to the Bernoulli-Laplace side of the 
spectrum. Either way, we hope in the future to work towards a more complete 
characterization of the mixing rates for various lattice path transposition chains. 

It is also not clear what consequences the new bound for Dyck paths might 
have for mixing on other Catalan structures. Many well-known Catalan struc¬ 
tures suggest natural Markov chains, but it is striking how different they can 

2 One example is the lattice path matroid C[A,B\ where A = (—b ) 71 and B = (H—) n . 

which is equivalent to the random walk on the n-cube slowed down by a factor of n. 
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be—even though there are natural bijections between classes of Catalan struc¬ 
tures, they frequently do not preserve natural notions of distance, making stan¬ 
dard comparison arguments tricky. There seems to be no obvious way to lever¬ 
age even a hypothetical O(nlogn) mixing time result for Dyck transpositions 
to yield improved bounds for any other Catalan chain. It is worth noting that 
the classical (2 n + 1-to-one many) Chung-Feller cyclic lemma suggests project¬ 
ing the chain consisting of all transpositions on the (Bernoulli-Laplace) space of 
{±l} 2ra+1 , with (n + 1) l’s and n -l’s, onto the set of Dyck paths. This projec¬ 
tion chain is in fact a Markov chain, as was observed by |3J , and as such inherits 
a lower bound on the spectral gap of order 1/n, from that of the Bernoulli- 
Laplace chain; however, many of the moves in the resulting projection chain are 
not particularly natural when viewed from the point of Dyck paths, and once 
again comparisons do not seem to help very much. 

An alternate approach would be to come up with a general scheme to bound 
mixing times which could be modified for any Catalan structure. For example, 
it is well-known that the Catalan numbers satisfy the recurrence 

n 

i= 1 

and this is evidenced by a recursive structure in almost every combinatorial 
realization of the Catalan sequence. It is common to exploit recursive structure 
in bounding mixing times, and indeed both the triangulation bound of m and 
the balanced matroid bounds of mil uni exploit different recursive structures 
of those objects. Unlike the recurrence above, the recursive structures used in 
those papers seem very particular to the Catalan realizations at hand, and do 
not generalize well to other Catalan chains. 

Finally, we have mentioned above the conjecture that the basis exchange 
walk is rapidly mixing for all matroids. As evidence, it would be interesting to 
show rapid mixing for some class of matroids which is not a subclass of balanced 
matroids. One promising candidate might be the class of transversal matroids, 
for which the basis exchange walk seems closely related to the well-studied 
problem of walks on bipartite graph matchings. 
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